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Abstract 

We perform a canonical, reduced phase space quantisation of General Relativity by Loop Quantum 
Gravity (LQG) methods. 

The explicit construction of the reduced phase space is made possible by the combination of 
1. the Brown - Kuchar mechanism in the presence of pressure free dust fields which allows to 
deparametrise the theory and 2. Rovelli's relational formalism in the extended version developed by 
Dittrich to construct the algebra of gauge invariant observables. 

Since the resulting algebra of observables is very simple, one can quantise it using the methods 
of LQG. Basically, the kinematical Hilbert space of non reduced LQG now becomes a physical 
Hilbert space and the kinematical results of LQG such as discreteness of spectra of geometrical 
operators now have physical meaning. The constraints have disappeared, however, the dynamics 
of the observables is driven by a physical Hamiltonian which is related to the Hamiltonian of the 
standard model (without dust) and which we quantise in this paper. 



1 Introduction 

The objects of ultimate interest in a field theory with gauge symmetry are the gauge invariant 
observables. There are two major approaches to the canonical quantisation of such theories. In 
the so called Dirac approach one first constructs Hilbert space representations of gauge variant non 

*gieskri@aei. mpg.de 

'l'thiemann@aei.mpg.de,tthiemann@perimeterinstitute.ca 



1 



observables and then imposes the vanishing of the quantised version of the classical gauge symmetry 
generators (constraints) as a selection principle for physical states. The associated physical Hilbert 
space then hopefully (if there are no anomalies) carries a representation of the observable algebra. 
In the so called reduced phase space approach one first constructs the classical observables and then 
directly looks for representations of that algebra. 

The advantage of the Dirac apporoach is that the unreduced phase space of non observables is 
typically a smooth (Banach) manifold so that the algebra of non - observables is sufficiently simple 
and representations thereof are easy to construct. Its disadvantage is that one has to deal with 
spurious degrees of freedom which is the possible source of ambiguities and anomalies in the gauge 
symmetry algebra. The advantage of the reduced phase space approach is that one never has to 
care about kinematical Hilbert space representations. However, its disadvantage is that the reduced 
phase space typically no longer is a smooth manifold turning the induced algebra of observables so 
difficult that representations thereof are hard to find. 

The reduced phase space of General Relativity with standard matter is hard to construct explicitly. 
However, on can combine two independent recent developments in order to make progress: 

On the one hand. Brown & Kuchaf have shown in a seminal paper [1] that there is hope to 
construct observables if one adds pressure free dust to the theory. This is because one can then write 
the constraints in deparametrised form0. 

On the other hand, there is Rovelli's relational formalism [2] for constructing observables which 
we need in the extended form developed by Dittrich p]. With this formalism one can write the 
observables as an infinite series Fj-t in terms of powers of so called clock variables T and with 
coefficients involving multiple Poisson brackets between constraints C and non observables / such 
that the series is (formally) gauge invarianlH. Remarkably [3l H] , the map Ft : f ^ Fj-t is a Poisson 
homomorphism between the algebra of non observables / and the algebra of observables with respect 
to a certain Dirac bracket (which is uniquely determined by the constraints and the functions T). 

Now usually Dirac brackets make the Poisson structure so complicated that one cannot find 
representations thereof. However, as observed in [1], if the system deparametrises, if one uses as 
clocks T the configuration variables conjugate to the momenta P m. C = P + H and if one considers 
functions / which do not depend on T, iijthen Ft becomes a Poisson bracket isomorphism. Moreover, 
the functions H m C = P -\- H become physical, conserved Hamiltonian densities which drive the 
physical evolution of the observables. This implies that a reduced phase space quantisation strategy 
becomes available, since to find representations of the Fj^t is as easy as for the /. The only non 
trivial problem left is to find representations which support the physical HamiltoniaiJ^. 

In [5] these two independent observations were combined and the algebra of classical physical 
Observables was constructed explicitly by adding a general scalar field Lagrangian without potential 
to the Einstein - Hilbert and standard model Lagrangian. It turns out that among the, in principle, 
infinite number of physical Observables there is a unique, positive Hamiltonian selected. 

In P, [7] that framework was further improved by using as specific scalar field the pressure free 
dust of Brown & Kuchaf. The corresponding Hamiltonian is positive, reduces to the ADM energy far 

"'^ Given a system of constraints Ci on a phase space, deparametrisation means that one can find local coordinates in 
the form of two mutually commuting sets of canonical pairs {q°',Pa), {T^ , ttj) such that the constraints can be written 
in the locally equivalent form Ci = i:i -\- Hi where the Hj only depend on the (q'^^pa)- 

^It is manifestly gauge invariant in an open neighbourhood of the phase space if the series converges with non zero 
convergence radius which has to be checked. 

•^This is no loss of generality because P can be eleminated in terms of the other degrees of freedom via the constraints 
and T is pure gauge. 

*The caveat is that the deparametrisation and thus the reduced phase space quantisation is generically only locally 
valid in phase space. Thus, the globally valid Dirac quantisation programme should be developed further in parallel. 
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away from the sources and to the standard model Hamiltonian on flat space. It generates equations of 
motion for the observables associated to the non dust variables that are in agreement with the Einstein 
equations for the system without dust, up to small corrections which originate from the presence 
of the dust. In particular one can develop a manifestly gauge invariant cosmological perturbation 
theory to all orders which was shown to reproduce the linear order as developed by Mukhanov, 
Feldmann and Brandenberger |[8|. The dust serves as a material reference system which we couple 
dynamically as fields rather than assuming the usual test observers in order to give the Einstein 
equations (modulo gauge freedom) the interpretation of evolution equations of observable quantities. 
This leads to in principle observable deviations from the standard formalism which however decay 
during the cosmological evolution. 

In this paper we quantise the algebra of observables constructed in [6]. Actually there is not 
much to do because that algebra is isomorphic to the Poisson algebra of General Relativity plus 
the standard model on M x 5 where S is the dust space manifold. Hence we can take over the 
kinematical Hilbert space representation that is used in Loop Quantum Gravity (LQG) [9|,[Tn]. For 
recent reviews on LQG see [H], for books see [12]. One may object that this representation is less 
natural here than in usual Dirac quantised LQG where it is uniquely selected on physical grounds 
p!3l [H] , namely one wants to have a unitary representation of the spatial diffeomorphism group of 
the coordinate manifold A" which is a gauge group (passive diffeomorphisms) there. Since all our 
observables are gauge invariant, we have no diffeomorphism gauge group any longer, hence that 
physical selection criterion is absent. However, it is replaced by a different one: It turns out that the 
physical Hamiltonian has the diffeomorphism group of the dust label space as symmetry group. These 
diffeomorphisms change our observables, they are active diffeomorphisms since they map between 
physically distinguishable dust space labels. Thus we may apply the same selection criterion. 

Now the interesting remaining question is whether that representation allows us to define the 
quantised version of the physical Hamiltonian. Maybe not surprisingly, it turns out that the same 
techniques that allowed to construct the quantum Hamiltonian constraint [TS] and the master con- 
straint [16] in usual Dirac quantised LQG can be used to define the quantised physical Hamiltonian. 
This operator is positive, hence symmetric and upon taking its natural Friedrich extension, it be- 
comes self - adjoint. In order to preserve its classical, active, spatial diffeomorphism symmetry it 
turns out that one has to define it in such a way that it preserves the graph of a spin network 
function that it acts on. The techniques developed in can now be applied to show, using the 
semiclassical states introduced in [18], that the physical Hamiltonian has the correct semiclassical 
limit on sufficinently fine graphs. In fact, in order to get rid of the graph dependence one can use the 
generalisation of LQG to Algebraic Quantum Gravity [17j. This casts quantum gravity completely 
into the framework of (Hamiltonian) lattice gauge theory [191 ISQ] with one crucial difference: There 
is no continuum limit to be taken because we are in a background independent theory with active 
diffeomorphisms as symmetries. 

The attractive feature of this reduced phase space approach is that we no longer need to deal 
with the constraints: No anomalies can arise, no master constraint needs to be constructed, no 
physical Hilbert space needs to be derived by complicated group averaging techniques. We map 
a conceptually complicated gauge system to the conceptually safe realm of an ordinary dynamical 
Hamiltonian system. The kinematical results of LQG such as discreteness of spectra of geometric 
operators now become physical predictions. This is a concrete implementation of the programme 
outlined for the full theory in and generalises the reduced phase space techniques recently adopted 
for the Loop Quantum Cosmology (LQG) truncation of LQG [22], [231 1211 [25] which is a toy model 
for the cosmological sector of LQG, to the full theory. 

It "remains" to analyse the physical Hamiltonian in detail since it encodes the complete dynamics 
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of General Relativity coupled to the standard model. The following tasks should be addressed in the 
future: 

1. Vacuum and spectral gap 

For a start we notice that the physical Hamiltonian does not depend explicitly on an external 
time parameter. Our Hamiltonian system which dynamically couples geometry and matter is 
a conservative system. This is in contrast to QFT on curved and in particular time dependent 
background spacetime metrics where one quantises matter propagating on an externally given 
background geometry. The Hamiltonian of that QFT is not preserved and thus even the notion 
of a ground state or vacuum as a lowest energy eigenstate becomes time dependent which leads 
to constant particle creation problems etc. ^26j. In our approach the notion of a vaccum state 
would not suffer from those problems. This appears as a conceptual improvement although of 
course the lowest eigenvalue of the Hamiltonian could be vastly degenerate. Also, the minimum 
of the spectrum of the Hamiltonian might not lie in its discrete (more precisely, pure point) 
part so that the "ground state (s)" would not be normlaisable. 

2. Scattering theory 

With a physical Hamiltonian H at our disposal we can in principle perform scattering theory, 
that is, we can compute matrix elements of the time evolution operator U{t) = exp(2rH). 
The analytical evaluation of those matrix elements is of course too difficult but as in ordinary 
QFT we may use Fermi's golden rule and expand, for short time intervals r, the exponential 
as U{t) = 1-)-/^ + ir H +0(r)^. The matrix elements of H seem hopeless to compute because it 
involves square roots of a positive self adjoint operator for whose precise evluation we would 
need the associated projection valued measure which of course we do not have. However, since 
in scattering theory initial and final states are excitations over a ground state which we do not 
know exactly but presumably can approximate by kinematical coherent states, one can invoke 
the technique developed in to expand the square root of the operator around the square 
root of its exectation value. We will do this in a future project. Of course there are issues to 
be resolved such as those of the existence of asymptotic states [27] and how one implements 
them in our formalism, see e.g. [2^ for some basic ideas. 

3. Anomalies 

As already mentioned, the Hamiltonian H has a huge symmetry group of which Diff(iS) is a 
subgroup and it is easy to implement this symmetry at the quantum level. However, there is 
another infinite classical, Abelian symmetry group JV which is generated by the Hamiltonian 
density functions H{a) and in terms of which the Hamiltonian reads H = d^a H{a). 
Classically one has {H{a),H{a')} = which of course implies classically that {H{a), H}=0. 
The Lie algebra of the total classical symmetry group thus consists of infinitesimal active 
diffeomorphisms and infinitesimal transformations generated by the H{a). The latter form an 
Abelian Poisson ideal and thus J\f is an Abelian invariant subgroup in the total symmetry group 
which hence is a semidirect product C5 = A/" xi DiS{a). Presumably, in the naive quantisation 
of H that we consider as a preliminary proposal in this paper, the latter symmetry is explicitly 
broken, or anomalous although semiclassically it is preserved. In order to reinstall it, one can 
try to make use of renormalisation group techniques associated to so called improved or perfect 
actions [29] . 

4. Lattice numerical methods 

It transpires that within the framework proposed here many of the conceptual problems of 
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canonical quantum gravity have been solved and the technical tasks have been simplified and 
reduced to a detailed analysis of the operator H, of course, at the price to have introduced 
additional, albeit unobservable, matter as a material reference system and a possibly only 
locally (in phase space) description. Since H is a complicated operator which is is formulated 
in terms of lattice like variables especially in the AQG version, it is natural to use Monte Carlo 
methods in order to study the operator numerically. 

5. QFT on curved spacetimes and standard model 

It is widely accepted that the framework of QFT on curved spacetimes |26] should be an 
excellent approximation to quantum gravity whenever the metric fiuctuations are small. In 
particular, when the background spacetime is Minkowski, then the standard model must be 
reproduced. Besides that, one would like to see whether our background independent lattice 
theory which is manifestly UV finite and non perturbative can explore the non perturbative 
sector of the standard model such as QCD. Another interesting question is whether our explic- 
itly geometry - mater coupled system can lead to an improved understanding of the Hawking 
effect due to the possibility to take care of backreaction effects. 

6. Effective action, universality, ambiguities 

Our framework presents a canonical quantisation of the field theory underlying the Einstein 
Hilbert action plus standard model action. Now computations within perturbative QFT and 
also string theory suggest that the effective actioijfl for gravity is an extension of the Einstein 
- Hilbert Lagrangian by higher derivative terms and an often asked question is whether one 
should not quantise these more general actions. There are several remarks in order: 

A. The effective action is a complicated, often even non local, action which takes care of all 
higher loop diagrammes obtained from a simple bare action. It looks like a classical action 
but it actually encodes all quantum fluctuations. Therefore it is inappropriate to quantise 
that classical action anew, it would not produce the same quantum theory as the bare 
action. 

B. Still one could just add all possible higher derivative terms from the outset. While one can 
canonically quantise such theoires by the Ostrogradsky formalism, this leads in general 
to a drastic increase in the number of degrees of freedom [31] due to the appearance of 
higher time derivatives. 

C. In the Euclidean formulation of QFT on Minkowski as a path integral one entertains a 
related (Wilson) notion of effective action as the action that one obtains when integrating 
out degrees of freedom labelled by (in Fourier space) momenta above a certain energy 
sea 1^. This also produces various higher derivative terms at lower energies as compared 

^There are several loosely equivalent definitions for the effective action. The notion we mean here is the following: 
Consider first a renormalisable theory. Given a defining action with a finite number of finite but unknown couplings 
and masses (parameters) one can perform perturbation theory and discovers, within a given regularisation scheme, 
that the parameters are to be altered by functions of the distance cutoff which diverge in the limit of vanishing cutoff 
in order to avoid singularities in loop diagrammes. If one does this order by order then one ends up with the so called 
bare action which produces finite higher loop diagrammes to all orders. The effective action is a vehicle that produces 
the same scattering amplitudes or n— point functions as the bare action but of which one only needs to compute tree 
diagrammes (no loops). The definition for a non renormalisable theory such as gravity is the same, just that then 
number of parameters is infinite. In renormalisable theories a finite number of experiments is sufficient to fix the 
unknown parameters while non renormalisable theories have no predictive power. 

^That energy scale has nothing to do with a perturbative cutoff, we are talking here about an already well defined 
theory. 
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to the bare action which is defined at infinite energy. Now the couphngs of the bare action 
also also are in principle unknown, however, for many theories that does not matter due 
to a phenomenon called universality: The couplings of the higher derivative terms depend 
on the energy scale and a coupling is called relevant, marginal or irrelevant respectively if 
it grows, remains constant or decreases in the low energy limit. A universal theory is such 
that all but a finite number of the couplings are irrelevant. One may ask whether one 
can see universality also in the canonical formalism, however, there are several obstacles 
in answering this question. First of all, the Euclidean formulation uses a Wick rotation 
which is only possible for background dependent theories where the background has a 
presentation with an analytic dependence on the time coordinate. In quantum gravity the 
metric becomes an operator, hence Wick rotation and therefore a Euclidean formulation 
is not possible. One should therefore define the Wilsonian effective action directly in the 
Hamiltonian (Lorentzian) formulation, however, that has not been done so far. 

It seems to us that in order to make progress on this kind of questions one should first try 
to define a Hamiltonian notion of effective action, see [33j for a possible direction. Then, if 
the symmetry arguments mentioned under [3.] are insufficient in order to fix the quantisation 
(discretisation) ambiguities in the definition of H, possibly universality studies may lead to 
further understanding. 

7. Singularity avoidance 

In quantum gravity we expect or want to resolve two types of singularities: First, QFT kind of 
short distance singularities which come from the fact that in interacting field theories one has to 
deal with products of operator valued distributions. Secondly, classical General Relativity kind 
of singularities which are simply a feature of the Einstein equations to predict that generically 
spacetimes are geodesically incomplete. An analytical measure for such spacetime singularities 
are typically divergences of curvature invariants. 

Now as shown in [15], UV type of singularities are absent at the non gauge invariant level, 
specifically, the quantum constarints are densely defined. In [21] it was discovered, in the con- 
text of usual LQG that expectation values of non gauge invariant curvature operators with 
respect to non gauge invariant coherent states that are peaked on a classically singular (FRW) 
trajectory remain finite as one reaches the singularity, thus backing up the much more spec- 
tacular results of [221 ISSl [211 [25] which are at the level of the physical Hilbert space albeit for 
a toy model and not the full theory. 

While these are encouraging results, they are at the kinematical level only and thus are in- 
conclusive. However, with the technology developed in this paper we can transfer both results 
literally and with absolutely no changes to the physical Hilbert space. As far as the spacetime 
singularity resolution is concerned, this is still not enough because the coherent states that 
we are using, while being now physical coherent states, they are not adapted to the physical 
Hamiltonian and thus may spread out under the quantum dynamics generated by U{t). In 
other words, given gauge invariant initial data m(0) and a coherent state ipQ that we prepare at 
r = and which is peaked on m(0), it may be that after short time r the state U{t)iPq is very 
different form the state ipr which is peaked on the classical trajectory r i— *• m{T). Therefore, 
in order to come to conclusions one should rather study expectation values with respect to the 
states U{t)%1)q rather than ip^-. In addition, one should try to construct dynamical coherent 
states for which such a spread does not happen. However, this is a difficult task already for 
the anharmonic oscillator. 
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The plan of the paper is as follows: 

In section two we review the essentials of O E] in order to make this article self - contained. This 
will lead to the reduced phase space and the classical physical Hamiltonian. 

In section three we quantise the reduced phase space using methods from LQG and obtain the 
physical Hilbert space almost for free. Then we implement the physical Hamiltonian on that Hilbert 
space. We do this both for LQG and the AQG extension. 

In section four we summarise and conclude. 



2 Review of the Brown — Kuchaf and relational framework 
2.1 Brown — Kuchar Lagrangian 

In [1] Brown and Kuchar add the following Lagrangian to the Einstein - Hilbert and standard model 
Lagrangian on the spacetime manifold M 



Sd = -1 [ d'X v^d^ p [g'^'' U,U, + 1] (2.1) 



where the one form U is defined hj U = —dT + WjdS^ and the index j takes values 1, 2, 3 while /x, v 
take values 0, 1, 2, 3. The action is a functional of the fields p, g^^, T, , Wj. Here T, have 
dimension of length, Wj is dimensionless and thus p has dimension cm~^. 

As shown in [HIE], in performing the Legendre transformation of fl2.ll) according to the 3+1 split 
of M = R X into time and space one introduces momenta P, Pj, I, P conjugate to T, S\ p, Wj 
respectively next to the momenta P°*, p, pa conjugate to qat, n, n°- respectively one encounters 
several primary constraints. Here one has introduced a foliation of M, that is, a one parameter 
family of embeddings t \^ Xt : X ^ Xt where Xt are the leaves of the folitation and the coordinates 
on X are denoted by x", a = 1,2, 3. The vector field dtX^ = nn^ + n°-X^^ can be dexompsed in 
components normal and tangential to the laves where is the future oriented normal. The functions 
n, n"' are the usual lapse and shift functions and qab = QtivX^^X^ defines the three metric intrinsic to 
X . The afore mentioned primary constraints are 

Z =: / = 0, := P = 0, Zj := Pj + PWj = 0, z := p = 0, Za := Pa = (2.2) 

The stability analysis of these constraints with respect to the corresponding primary Hamiltonian 
leads to the following secondary constraints 



P 



2 



c^ot ^ c + c^, = - [ + p^d^(l + q^'UM)] 

^ PVdet(g) 

cf = c, + cf, c^ = P[T,-W,S^J 

c = , C-^ + Vd^){l + q'''UaU,)] (2.3) 

and six more equations which can be solved for the Lagrange multipliers corresponding to constraints 
Z^ , Zj and which we do not display here. Here Ua = —Ta + WjS^^ = —c^/P and c, cf respec- 
tively are the contributions of geometry and standard matter to the usual Hamiltonian and spatial 
diffeomorphism constraint respectively. 



7 



The stability analysis of the secondary constraints with respect to the primary Hamiltonian 
which is a linear combination of the constraints fl2.2p and the first two constraints in (12.31) reveals 
that there are no tertiary constraints. Moreover, the classification of the sets of constraints into first 
and second class shows that the constraints 2;, Za, c*°*, c*°* are first class while, roughly speaking, 
the pairs {Z,c), {Zj, Z^) form second class constraints with non degenerate matrix formed by their 
mutual Poisson brackets. Hence, to proceed, one passes to the corresponding Dirac bracket and 
solves the second class constraints explicitly by setting 

/ := 0, := 0, W, := := -^l^[q-^UaU, + 1] (2.4) 

^ Vdet(g) 

Fortunately, the Dirac bracket reduced to the geometry variables g^b, p""^ and the remaining matter 
variables is identical to the original Poisson bracket. 

After using (12.41) and solving z = Zahy identifying lapse and shift as Lagrange multiplicator 
functions respectively we are left with the first class constraints 



Cf = Ca + Cf , = PTa + P,S^,a (2-5) 

In principle we could have chosen the other sign to solve the quadratic equation for p in (12.41) but the 
detailed analysis in [6] reveals that the other choice would produce the Einstein equations with the 
wrong sign in the limit of vanishing dust fields. In particular one must choose p, P < so that the 
additional matter enters with negative sign into the Hamiltonain constraint. This has the important 
consequence that c > thus enables close to flat space solutions. 

As far as the physical interpretation of the additional matter is concerned we just mention that 
its Euler Lagrange equations imply that the vector field = g^^Up is a geodesic in affine parametri- 
sation, that the fields Wj, are constant along the geodesic and that the field T defines proper time 
along each geodesic. It follows that = =const. labels a geodesic while T = r =const. is an 
affine parameter along the geodesic. Furthermore, its energy momentum tensor is that of a perfect 
fluid with vanishing pressure and negative energy densitju, hence it is pressure free phantom dust. 
It serves as a dynamical, material reference system which also plays the role of a phantom in the 
literal sense because it is not directly visible in the final picture while leaving its fingerprint on the 
dynamics. 



2.2 Brown — Kuchar Mechanism 

The observation of Brown and Kuchar was that the constraints (12. 5 p can be written in deparametrised 
form. This holds in more general circumstances, namely whenever we consider scalar fields without 
potential and mass terms as pointed out in [5]. The observation consists in the fact that the only 
appearance of T, in c*°* is in the form . However, this means that using c*°* = we may write 
(12. 5p in the equivalent form 

c*"* = c + c^, = -^/P^Tq^^b 

C = ca + cf , = PT, + P,S% (2.6) 

^We are not violating any energy conditions because we still require that the energy momentum tensor of observable 
(standard) matter plus dust satisfies the energy conditions. In fact, it would be sufficient if the energy conditions are 
satisfied by the standard matter alone because in the final analysis the dust completely disappears while the equations 
of motion for observable matter and geometry assume their standard form plus small corrections, see [Bl|7]- 
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where equivalent means that fl2.5p and (12.61) define the same constraint surface and the same gauge 
invariant functions. 

We can now solve the first equation in (12. 6p for P, remembering that P < and the second 
equation for Pj, making the assumption that the matrix S^^ is everywhere non degenerat^ with 
inverse Sj. The result is 

cf = p^ + h„ h, = s;;[ca-hTa] (2.7) 

In solving (12. 6p in terms of P we find at an intermediate step that = — q°'''CaCb. Hence, while 
the argument of the square root in (12. 7p is not manifestly positive, it is constrained to be positive. 

Notice that the function h is independent of S^,T while hj still depends on both. Hence, we 
have achieved only partial deparametrisation. However, this will be sufficient for our purposes. An 
important consequence is that the constraints in the form (12. 7p are mutually Poisson commuting. 
This follows immediately from an abstract argument^ [34] , although one can also verify this by direct 
computation JT\. This implies in particular that the h{x) are mutually Poisson commuting while the 
h{x) do not Poisson commute with the hj{y) and neither do the hj{y) among each other. 

2.3 Relational framework 
2.3.1 General theory 

We first consider a general system with first class constraints Cj with arbitrary index set X and later 
specialise to our situation. 

Consider any set of functions on phase space such that the matrix defined by the Poisson bracket 
entries M/ := {C/,T"^} is invertible. Consider the equivalent set of constraints 

c'r.= Y.lM-YjCj (2.8) 
J 

such that {Cj, Tj} ^ 5/ where ~ means = modulo terms that vanish on the constraint surface. Let 
Xj be the Hamiltonian vector field of Cj and set for any set of real numbers 

Xp:=J2P'Xj (2.9) 

For any function / on phase space we set 

oo ^ 

a^if) ■■= exp(X^) ■f = Y.^_^^p-f (2.10) 

n=0 

Now let be another set of real numbers and define 

0/(r) := [aMUiT)=r (2.11) 
^This is a classical restriction of the same kind as det{q) > 0. 

^ The constraints (|2.7p are first class. Hence their Poisson brackets are linear combinations of constraints. Since 
the constraints are linear in the momenta P, Pj , their Poisson brackets are independent of P, Pj . Therefore we can 
evaluate the linear combination of the constraints that appear in the Poisson bracket computation in particular at 
P = -h, Pj = -hj. 
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where ^/^(T) = r means apiT^) = for all I. As one can check, ap{T^) ^ + (3^ so that fl2.1ip is 
weakly (i.e. on the constraint surface) equivalent to 

Of{T) := Mf)h=,-T (2.12) 

Notice that after equating /3 with t — T, the previously phase space independent quantities P become 

phase space dependent, therefore it is important in fl2.12p to first compute the action of with /3 
treated as phase space independent and only then to set it equal to r — T. 
The significance of (12.121) lies in the following facts: 

1. The functions Oj(r) are weak Dirac observables with respect to the Cj, that is 

{Ci,Of{T)}^0 (2.13) 

This remarkable property is due to the key observation that the Xj weakly commute [3l H] . 

2. The multi parameter family of maps : / 0/(r) is a homomorphism from the commutative 
algebra of functions on phase space to the commutative algebra of weak Dirac observables, both 
with pointwise multiplication, that is 

Ofir) + Ofir) = Of^fir), O^r) Of^r) ^ Ojfir) (2.14) 

The linear relation is obvious, the multiplicative one follows from the fact that 

Mff) = ■ //' = e""^ . ffe-""^ ■ 1 = [e^'' ■ /e-^^][e^''/'e-^''] (2.15) 

where we used the identity 

[e^../e-^.] = ^ _[X,,/](„) (2.16) 

n=0 

and where Xp, f respectively are considered as derivation and multiplication operators respec- 
tively on the algebra of functions on phase space so that [X/j, /] = X/3 ■ /. Here [X, /](o) = 

/, [X,f]{n+1) = [X, [X, /](„)]. 

3. The multi parameter family of maps : f ^ Of{T) is in fact a Poisson homomorphism with 
respect to the Dirac bracket {., .}* defined by the second class system Cj, T"^, that is 

{Ofir),Ofir)} ^ {Ofir),Of>ir)r^O{fj>yir) (2.17) 

where the Dirac bracket is explicitly given by 

{/, f'V = {f, n - {f, Cj}[M-'Yj{T', /'} + {/', Cr}[M-'Yj{T', /} (2.18) 

Here we have used in the first step that both Of{T), Oflr) have weakly vanishing brackets 
with the constraints. Relation (12.171) follows from the fact that the map a^j is a Poisson 
automorphism on the algebra of functions on phase space and the Poisson bracket must be 
replaced by the Dirac bracket because in evaluating {0/(r), 0//(r)} we must take care of the 
fact that P = T — T is phase space dependent. See |1] for the explicit proof. 

The interpretation of Of{T) is that it is a relational observable, namely it is the value of / in the 
gauge P = T — T. 
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2.3.2 Specialisation to deparametrised theories 

For deparametrised theories it is possible to find canonical coordinates consisting of two sets of 
canonical pairs (P^,T/) and {q"',Pa) respectively (where the Poisson brackets between elements of 
the first and second set set vanish) such that the constraints Cj can be rewritten in the equivalent 
form 

Ci = Pi + hiiq'',Pa) (2.19) 

that is, they no longer depend on the variables T^. This is a very special case and most gauge 
systems cannot be written in this form. Even with dust General Relativity is a priori not of that 
form, however, we will reduce it to that form with an additional manipulation below. 
The simplifications that occur are now the following: 

A. We obviosly have 

M/ = {C/,T*^} = (5/ (2.20) 
therefore Cj = Cj and we do not have to invert a complicated matrix. 

B. By the same argument as in the footnote after (12. 7p we have {C/, Cj} = identically on the 
full phase space, not only on the constraint surface which of course implies that = 0, 
the Hamiltonian vector fields of the constraints are mutually commuting. It also follows that 
{hj^hj} = and thus {Ci,hj} = for all I, J which means that the hj are already Dirac 
observables. 

These simplifications mean that all the previous weak equalities become strong ones, i.e. identities 
on the full phase space. The Dirac observable associated to Tj 

OtAt) = MT%,iT)=r = r' (2.21) 

is simply the constant (on phase space) function r^. The momenta Pj are already Dirac observ- 
ables, however they can be expressed in terms of q°',Pa via the constraints. Moreover, since is a 
homomorphism we have on the constraint surface 

Pi = Op,{t) = -OhAr) = -/^/(0,a(r), Op.(r)) =: -Hj (2.22) 

In fact we have hj = Hj because hj is already a Dirac observable. 

The reduced phase space (where the constraints hold and where the gauge transformations have 
been factored out) is therefore coordinatised by the functions 

Q%T) = OAr), Pa{r) = 0,Sr) (2.23) 

and in what follows we concentrate on functions / which only depend on q°',Pa- On such functions 
the Dirac bracket reduces to the Poisson bracket since {T/, /} = for all I. Therefore the reduced 
map O'^ : / i—> 0/(r) is now a multi - parameter Poisson automorphism with respect to the Poisson 
bracket. In particular we note 

{PaiT),Q\r)} = {OpAr),OAr)} = 0{,„,,.}(r) = 0,.(r) = 6^ (2.24) 

which means that the reduced phase space has a very simple symplectic structure in terms of the 
coordinates Pa '■= -Pa(O), Q'^ := <5"(0) which in fact form a conjugate pair. It is this fact which 
makes reduced phase space quantisation feasible as observed in [1]. 
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It seems that we have triviahsed everything. However, this is not the case as we must interprete 
the r dependence of our observables. We notice first of all that on functions / independent of , Pj 
formula fl2.12p reads explicitly 

Ofir) = arU) = exp(X.) " / = E " ^ (2-25) 

n=0 

where Xr is the Hamiltonian vector field of the function Hr = {t^ — T^)Hj. Here we have used that 
the Xj on / reduce to the Hamiltonian vector field of hj and since hj is independent of Pj we may 
write Of{T) in the above compact form. It is now a simple exercise to verify that [1] 

dOf{T)/dT' = {Hi,Of{r)} (2.26) 

which means that the strongly Abelian group of Poisson bracket automorphisms ar is generated by 
the "Hamiltonians" Hj. Thus, if we interpret the Tj as clocks then we have a multi - fingered time 
evolution with Hamiltonians Hj. 

In quantum theory then one would like to select a suitable one parameter family by prescribing 
functions t\s) in terms of a single parameter such that the associated Hamiltonian is positive and 
has preferred physical properties. 

2.4 The reduced phase space of General Relativity with dust 

Now we specialise to our situation which is a special case of the general theory. This has been previ- 
ously done in detail, including proofs, in ^ and was also reviewed in j6|. Here we summarise those 
results. 

As previously mentioned, the Hamiltonian constraints in 02. 7p are in deparametrised form, how- 
ever, the spatial diffeomorphism constraints are not. However, the idea is to exploit the fact that 
the constraints (12. 7p are mutually Poisson commuting so that one can perform the reduction of the 
phase space in two steps: First we reduce with respect to the spatial diffeomorphism constraint and 
then with respect to the Hamiltonian constraint. More precisely, consider arbitrary functions 
on X and denote by the Hamiltonian vector field of the function 

cf := I d^a (3^'{x)c'°\x) (2.27) 

where we have defined c^q^ = c*°*. Then for arbitrary functions t^{x) = t{x), t^{x) := a^{x) on X 
the general formula reads 

Of{r) = [ap{f)U^^T)=r, apif) = expiXfs) ■ f (2.28) 

where T°(x) = T{x), T\x) = S^{x). We readily compute that a^(T^(a;)) = T'^{x) + [3^{x) so that 

Of{r) = K(/)]/3=.-T (2.29) 

Now since S^{x) Poisson commutes with c*°*(?/) we may rewrite fl2.29p in the form 

Of{T) = Ma^if))]^^,__g)]f,o=,.T (2.30) 

It turns out that one can compute the inner argument of fl2.30p rather explicitly with an immediate 
physical interpretation for judicious choices of the functions a^{x). Namely, for any scalar function / 
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built from of T, P, qabiP"'^ and the matter of the standard model one finds explicitly that for constant 
functions 

M/(^))]/?=.-5 = /(^)5»=. (2-31) 

In other words, whatever the value of x at which the function / is evaluated, fl2.3ip evaluates it at 
the point at which S^{x) assumes the value Since we have assumed that S^^ is everywhere 
invertible and thus defines a diffeomorphism between X and the range of which is the dust space 
S, the value is unique. Formula (12.311) is proved explicitly in [6] and will not be repeated here. 
Thus, f l2.3ip takes a simple form if we choose as / one of the following functions on S 

T:=T,P=j, q,u := qabS^i fP' := f (2.32) 

where 

J := det{dS/dx) (2.33) 

as well as 

a] := a^S"), e\ := ^y^, i^ai ■= '4'ai, i'ai ■= 'ipai, 0/ := 0/, := y (2.34) 

for connections a[, electric fields ej, fermions ipai, '4'ai and Higgs fields 0/ with conjugate momentum 
TT^ of the standard model where / labels a basis in the Lie algebra of the appropriate gauge group, 
see [l5j for the canonical formulation of the standard model coupled to gravity including appropriate 
background independent Hilbert space representations. 

It is clear that the evaluation of the functions (I2.32p and (I2.34p at x^ is nothing else than the 
pull back of the corresponding fields to S under the inverse of the diffeomorphism : X ^ S. We 
will denote the corresponding tensor fields on S as in (12.320 and (12.340 . Notice that while these are 
scalars on X they are tensor densities of the same weight on S as they havj^ on X. In [H E] it is 
shown that one can arrive at the spatially diffeomorphism invariant functions (I2.32p and (12.341) also 
by symplectic reduction with respect to the spatial diffeomorphism constraint which is an alternative 
proof of the fact that canonical pairs without tilde on X are mapped to canonical pairs on S. For 
instance 

{fP\a),~qU^')} = ^^Sij'^)Sia, a') (2.35) 

where k = IGTrGNewton- This also shows that it is sufficient to consider constant rather than 
arbitrary functions. 

Returning to (I2.30p we see that it remains to compute 

Of{T,a):=[apo{f{a))]f,o=r-T (2.36) 

where / is now an arbitrary function of the spatially diffeomorphism invariant functions (I2.32p and 
(I2.34p . Now we can use the simplified theory of section [2.3.2l because c*°* is written in deparametrised 

^°This statement sounds contradictory because of the following subtlety: We have e.g. the three quantities 
P(x), P{x) — P{x) / J{x), P{a) — P{xa). On X, P{x) is a scalar density while P{x) is a scalar. Pulling back 
P{x) to S = S{X) by the diffeomorphism a i-^ S^^{a) results in P{a). But pulling back P{x) back to S results in 
the 507716 quantity P{(j)- Since a diffeomorphism does not change the density weight, we would get the contradiction 
that P{(j) has both density weights zero and one on S. The resolution of the puzzle is that what determines the 
density weight of P{x) on X is its transformation behaviour under canonical transformations generated by the total 
spatial diffeomorphism constraint c^°* = + where c^, are the dust and non dust contributions respectively. 
After the reduction of c^°*, what determines the density weight of P(cr) on S is its transformation behaviour under 
([ca + PT^a]Sj/J){xa) — Cj{a) + P{a)Tj{a) and this shows that P{(t) has density weight one. 
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form, i.e. it does not involve T, any longer. Actually, formula fl2.36p would be awkward for non 
constant functions r because it depends on 



cf = / d^x (r - T){x) c*°*(x) (2.37) 



„tot / ^ ^\ ^tot 

JX 

which is expressed on the space X rather than dust space S. However, for constant r (12.371) is the 
integral of a density of weight one and can then be written in the form 

cf = f £a (r - f ) (a) [P + h] (a) (2.38) 

where 



s 



a] 



c((^) = j(a;a) 

cM) = "-^M (2.39) 

Notice that e.g. c is just the pull back of c and that one simply has to replace every tensor without 
tilde by their pulled back image with tilde. Thus constant r is uniquely selected by the requirement 
that c*°* is spatially diffeomorphism invariant. 
It follows now from section 12.3.21 that 



n=0 

and that 



oo „ 
0/(r,a) = ^-{iJ„/(a)}(„), = £ a [r - T]{a) h{a) 



(2.40) 



—Of{a,T) = {H,Of{a,T)}, H=/ d'a h{a) (2.41) 



Since the h{x) are mutually Poisson comuting it follows that also the h{a) are mutually Poisson 
commuting so that 

H{a, t) := a^o(/i((T))^o=,_r = H^r) =: H{a) (2.42) 

is independent of r and already a Dirac observable. 

Notice that the physical Hamiltonian H is positive. It enjoys the following symmetries: Since 
it is an integral over a density of weight one it is invariant under diffeomorphisms of S. Notice 
that 5 is a label space for geodesies and not a coordinate manifold, hence in contrast to the passive 
diffeomorphism group Diff(A'), the group Diff(5) are active diffeomorphisms. In particular, it follows 
that 

{H,9(a)} = (2.43) 

which also is a consequence of having chosen constant r, in which case the physical Hamiltonian has 
a maximal amount of symmetry. Had we not chosen constant r then the physical Hamiltonian would 
not be a Dirac observable. 
This also implies that 

{H,C,(a)} = 0, C(a,r) := a^o(g,(a))^o=._r =: C,{a) (2.44) 
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is actually independent of r, although cj ^ Cj. Notice that 

H{a) = ,Jc{a,ry - Q^'^{a,r)C,{a)Ck{a), C{cx, r) := apo{d{a))po=,_T (2.45) 
The second symmetry of H is of course that 

{H,H{a)} = (2.46) 
Let us write for some scalar and vector test functions /, respectively 

H{f):= I d^a f{(T) H{(r), C{u) := I d"^ a {a) C j{a) (2.47) 
Js Js 

then 

{C{u),C{u')} = -kC{[u,u']) 
{C{u),H{f)} = -^H{u[f]) 

{Hif),H{f')} = (2.48) 

which shows that the symmetry generators generate an honest Lie algebra g in contrast to the Dirac 
algebra underlying GR as was pointed out already in [1] and further examined in [55]. That Lie 
algebra has a subalgebra generated by the C{u) and an Abelian ideal generated by the H{f), hence 
it is not semisimple. The corresponding Lie group = A/'xDiff (5) is therefore the semidirect product 
of the Abelian invariant subgroup A/" to which the H{f ) exponentiate and the active diffeomorphism 
group of dust space. 

2.5 Physical interpretation and comparison with unreduced formaUsm 

The symmetry algebra g and the associated conservation laws play a crucial role in showing [6J that 
the equations of motion for the canonical pairs of true degrees of freedom 

(Q,fc, P^'); (4, Ej); (vj/,,, ^^j); ($,, tf) (2.49) 

which are the images of the canonical pairs 

(qjk, (Sj, e{)- {tPaj, ^aj); (2.50) 

under a^o(.)^o=^_T at r = assume the standard form that they have in General Relativity without 
dust |36j, with two important modifications: First, in usual General Relativity without dust the 
equations of motion generated by the canonical Hamiltonian h{n, n) = c{n) + c(n) which is a linear 
combination of the smeared Hamiltonian constraint c{n) = d^xnc and spatial diffeomorphism 
constraint c{n) = d^xTf-Ca-, involve arbitrary lapse and shift functions n, ra'^ on X which are in- 
dependent of phase space. However, in our formalism lapse and shift functions become dynamical 
function^ on iS, namely N = C /H and A^-' = —Q^^Ck/H. Secondly, without dust we still have 
constraints c = = while we have energy - momentum conservation laws H = e, Cj = —ej 
where e, ej are arbitrary functions on S independent of r. This turns dynamical lapse and shift into 
a function of Q^^^ej/e. The functions e, express the influence of the dust on the other variables 
and are the price to pay for having a manifestly gauge invariant formalism rather than assuming non 
dynamical test observers that turn geometry and matter into observable quantities. 

This concludes the classical analysis and the review of [B]. 

^^This is similar in spirit to [37] where one replaces lapse and shift test fields by hand by phase space dependent 
functions, carefully chosen (via Witten spinor techniques that enter the proof of the gravitational positive energy 
theorem) so that the resulting Hamiltonian is positive, at least on shell. 
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3 Reduced phase space quantisation of General Relativity 



3.1 Hilbert space representation 

Let us summarise the result of the previous section: By using the relational formalism we can 
explicitly compute the reduced phase space of General Relativity with dust. It is identical to the 
unreduced phase space without dust with proper identification of X with S and of the gauge invariant 
canonical pairs (12.491) with the gauge variant canonical pairs 

{qab, p"'); (4, e'j); {^^i, M; i<Pi, (3-1) 

of geometry and standard matter. The constraints have disappeared, they have been solved and re- 
duced. Instead of a linear combination of constraints on the gauge variant phase space coordinatised 
by (13. ip which generates gauge transformations, there is a physical Hamiltonian (I2.4ip which gener- 
ates physical time evolution on the gauge invariant phase coordinatised by ( I2.49p . From the classical 
point of view one should now simply solve those equations in physically interesting situations. In 
P, [7j we have done this in the context of cosmological perturbation theory [6l [7j which is written 
in manifestly gauge invariant form. This not only reproduces the standard results [8] but also will 
allow us to investigate higher order perturbation theory without running into problems with gauge 
invariance. 

In the quantum theory we are looking for representations of the Poisson *— algebra generated 
by (12.491) which supports a quantised version of the Hamiltonian H. The selection of appropriate 
representations will be guided by the symmetry group unveiled in the previous section. First of all, 
since we consider fermionic matter we are forced to work with tetrads rather four metrics. We use 
the second order formalism as displayed in [15] (that is, we write the Einstein Hilbert Lagrangian in 
terms of the spin conection of the tetrad which involves second order derivatives rather than using the 
first order Palatini formalism) in order to avoid torsion. This means that we formulate the geometry 
phase space in terms of su(2) connections and canonically conjugate fields {Aj, Ej) rather than in 
terms of the ADM variables Qjk, P^^ where / is an su(2) index. This casts the geometry sector of 
the phase space into a SU(2) Yang - Mills theory description. Th price to pay is that there is an 
additional Gauss constraint on the phase space (which has been reduced only with respect to the 
Hamiltonian and spatial diffeomorphism constraint) given by 

Gi := djE{ + euxAjE^j^ + fermion terms (3.2) 

just as for the matter Yang - Mills variables (we assume that the Cartan Killing metric is always 
^JK by appropriate normalisation of the Lie algebra basis). 

The gauge field language suggests to formulate the theory in terms of holonomies along one 
dimensional paths and electric fluxes through two dimensional surfaces, just as in unreduced LQG. 
There one has a uniqueness result p^[T^ which says that cyclic representations of the holonomy - flux 
algebra which implement a unitary representation of the spatial diffeomorphism gauge group Diff(A') 
are unique and are unitarily equivalent to the Ashtekar - Isham - Lewandowski representation [^[TU]. 
In our case we do not have a diffeomorphism gauge group but rather a diffeomorphism symmetry 
group Diff(A') of the physical Hamiltonian H. This is physical input enough to also insist on cyclic 
Diff(5) covariant representations and correspondingly we can copy the uniqueness result. 

Thus we simply choose the background independent and active diffeomorphism covariant Hilbert 
space representation of LQG used extensively in P3] and we ask whether that representation supports 
a quantum operator corresponding to H. 
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3.2 Subtleties with the Gauss constraints 



Before we analyse this question in detail, we should mention a subtlety: When one rewrites the 
geometry and standard matter contributions c, Ca to the total Hamiltonian and and spatial diffeo- 
morphism constraint in terms of the gauge theory variables, one can do this is G invariant form 
(where G is the compact gauge group underlying the corresponding Yang Mills theory) only by in- 
troducing terms proportional to the Gauss constraint, see e.g. [12]. For instance, the contribution to 
the spatial diffeomorphism constraint of a Yang Mills field on the unreduced phase space is given by 



constants of the corresponding Lie algebra. The function (13.31) really generates Yang Mills gauge 
transformations, however, it is itself of course not Yang - Mills gauge invariant due to the term 
proportional to the Gauss constraint 



Likewise, the geometry contribution c^*^" to c contains a term proportional to g^i° [I2] (however, 
does not). As far as the definition of the complete constraint surface is concerned, one can drop the 
various Gauss law contributions to c, Ca since we impose the Gauss laws independently anyway. This 
gives an equivalent set of constraints which is such that c, Ca are manifestly invariant under Yang - 
Mills type of Gauss transformations. However, now the algebra of the c*°*, c*°* only closes up to a 
term proportional to the various Gauss laws. 

The question is now whether this spoils the argument that the constraints in the form (12.71) 
are mutually Poisson commuting. In fact, we only can conclude that their Poisson brackets are 
proportional to c*°*, c*°* and the various while they must not depend on the dust momenta 
P, Pj. This means that their Poisson brackets are proportional to a Yang Mills gauge invariant linear 
combination of Gauss constraints. Hence, indeed the constraints c*°*, c*°* are Abelian only on the 
constraint surface of the Gauss constraints. 

This poses the question which consequences this has for the formalism developed in the previous 
section. First of all, all relations that we have written there remain valid modulo terms proportional 
to the Gauss constraints. Secondly, the physical Hamiltonian is manifestly Yang - Mills gauge 
invariant, manifestly Diff(A') invariant and invariant modulo the Gauss constraints under H . 

The strategy that we adopt is the following. In the presence of gauge fields we actually work 
with the non Gauss invariant contributions to the spatial diffeomorpphism constraints as in (13.31) 
and with the non Gauss invariant contribution to c^*^" such that algebra of Hamiltonian and spatial 
diffemorphism constraints closes without involvement of the Gauss constraints. This makes the 
analysis of the previous section go through without modifications at the price that the physical 
Hamiltonian is not Gauss invariant. When we quantise it turns out that one can actually solve the 
various Gauss constraints explicitly by Dirac constraint quantisation. That is, the Hilbert space 
can be projected to the Gauss invariant subspace which has an explicitly known orthonormal basis 
given by the Gauss invariant spin network functions (and their analog for the gauge group of the 
standard model). Therefore, on the Gauss invariant Hilbert space one can actually replace the 
by Cf^^^ because the correction term proportional to the Gauss constraint vanishes on the Gauss 
invariant Hilbert space (upon appropriate ordering of the Gauss constraint operator to the right so 
that no commutator terms arise). Thus Cj is replaced by its Gauss invariant analog and similarly 
one can replace C by its Gauss invariant analog so that H and H become manifestly Gauss invariant 
operators and H should have the symmetry group C5 as well. 



where fl-^ 




YM 



dael + euKaie^K 



(3.4) 
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An alternative route would be to also reduce the phase space with respect to the Gauss constraints, 
possibly using the framework of [38] and references therein. 

3.3 Quantum Hamiltonian 
3.3.1 Sign issues and strategy 

Before we go into details we must worry about yet another issue: As we have seen in the classical 
analysis, the expression = — Q^^CjCk is constrained to be non negative. Actually we have 
seen this only for — q^^CaCh but as we showed 

(^2 - Q^%Ct)icr) = ([c2 - q'''caCk]/J){x^) (3.5) 

and J > bys assumption (we have imposed J ^ everywhere, hence either J > everywhere or 
J < everywhere by continuity and we choose the first option). However, on the full, reduced phase 
space — Q^^CjCk maybe indefinite. In the quantum theory we therefore should derive, roughly 
speaking, a self adjoint operator (valued distribution) for H'^{<j) and restrict the spectral resolution 
of the Hilbert space to the positive spectrum part. This has to be done for every a. This maybe 
impossible because the corresponding spectral conditions could be incompatible. However, as already 
ponted out by Brown and Kuchaf [1], if we indeed manage to quantise -ff^(cr) in such a way that 
they are mutually commutinj^ then the corresponding spectral projections commute and the above 
requirement is consistent. Unfortunately, not only may it be hard to achieve commutativity of the 
operators corresponding to the various if^(cr), moreover it will be hard to compute the corresponding 
projection valued measures. 

Therefore, as a first step, in this article we adopt the following strategy: Classically, in the 
interesting part of the phase space we have — Q^^CjCk > 0. Therefore on this part of the phase 
space we have trivially — Q^^CjCk = — Q^''CjCk\- Hence on that part of the phase space we 
have the identity 



H = ^\C^-Q^^C,Ck\ = ^j\{\C^ - Q^^Cfik] + - Q^^C,Ck\) (3.6) 

The virtue of this rewriting is that both expressions, which are identical on the physically interesting 
piece of the phase space, can be extended to the full phase space without becoming imaginary. In the 
second version, the function actually vanishes on the unphysical part of the phase space. In either 
form, the square root now makes sense in the quantum theory because its argument is now a non 
negative expression. 

We remark that a discussion of similar sign issues and whether one should allow states in the 
quantum theory which violate the classical positivity of H^{ci) which is enforced by a constraint of 
the form — = and where is not manifestly positive while surely is, can be found 
for instance in [52] • There the authors argue that one should allow negative energy states because 
otherwise one would exclude the tunneling effects into the classically not allows regions which, as 
we know from quantum mechanical experiments, do happen. What happens mathematically is that 
in the operator constraint method (Dirac approach) one quantises both P and as self - adjoint 
operators on the kinematical Hilbert space and then solves the quantum constraint. The elements 

^^More precisely, one has to demand that the projection valued measures E„ for the H^{(y) mutually commute in 
order to avoid domain questions. Notice that the Poisson commutativity of the -ff (cr) implies the Poisson commutativity 
of the H^{(j) and vice versa. 
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of the corresponding physical Hilbert space may have support in the classically not allowed region 
of the configuration space (where they typically decay rather than oscillate) so that the expectation 
value of = becomes negative. This is possible only because the operator corresponding to P, 
while being a quantum Dirac observable, does not descend to a self adjoint operator on the physical 
Hilbert space. In a strict reduced phase space quantisation one would have to restrict the physical 
Hilbert space to states which have support only in the classically allowed region of the phase space 
and this may well be the physically correct procedure. However, for the moment, as we do not yet 
have sufficient control over the spectrum of H^, we comply with the conclusion of [39] and do not 
make any restriction on the physical Hilbert space. 

Thus, in this article we therefore propose to quantise the first version of (13. 6p which is a classically 
valid starting poinlQ. We then adopt a naive quantisation strategy and are able to construct a well 
defined Hamiltonian operator. That quantisation not necessarily has the property that the quantised 
versions of the H'^{(y) are mutually commuting and therefore the operator constructed in this paper 
should only be considerd as a preliminary step. However, that operator has the following three 
properties: It is manifestly Gauss invariant, manifestly Diff(iS) covariant and has the correct classical 
limit in the sense of expectation values and fiuctuations with respect to coherent states. However, 
it maybe anomalous with respect to the group M. In fact, the absence of that anomaly would be 
mathematically equivalent to showing that the Dirac algebra of General Relativity is implemented 
non anomalously. We stress, however, that the gauge symmetries of General Relativity have been 
exactly taken care of in the reduced phase space approach. We are talking here about a symmetry 
group and not a gauge group. To break a local gauge group is usually physically inacceptable 
especially in renormalisable theories where the corresponding Ward identities find their way into the 
renormalisation theorems. However, it may or may not be acceptable that a physical symmetry is 
(spontaneouly, explicitly ...) broken. For instance, the explicit breaking of the axial vector current 
Ward identity in QED, also called the ABJ anomaly, is experimentally verified. 

In lack of a physical justification for why the M symmetry should be broken, we view that 
potential anomaly as an indication that the quantisation of the present paper has to be improved. In 
fact, since we are effectively working with a background independent lattice gauge theory, it is useful 
to adopt strategies from lattice gauge theory in order to restore symmetries on the lattice that are 
broken in a naive quantisation. It turns out that in fortunate cases one can restore the symmetry 
by making the operator quasi non local. That is, in addition to next neighbour interactions one has 
to consider next to next neighbour interactions etc. which makes the action non local, however the 
coefficients of those additional interactions decay exponentially with the lattice distance. See e.g. 
|29j and references therein. 

We consider the completion of this step as a future research programme. In the course of that 
analysis we might even be able to fix the quantisation (discretisation) ambiguities, i.e. the coefficients 
in front of the various n-th neighbour contributions. 

With this cautionary remarks out of the way, we can now consider a naive quantisation of the 
Hamiltonian which is strongly guided by analogous techniques developed for the Hamiltonian and 
Master constraint of unreduced LQG [ISl [IS] so that these constructions are also helpful in the 
present reduced phase space approach. 

3.3.2 Naive Quantisation 



^•^A similar strategy was adopted for the quantisation of the volume in LQG: Classically we have det(q) — det(£') > 
but in order to give meaning to det(q) in the quantum theory we must start from yj\ det(i?)|. 
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3.3.2.1 Classical regularisation We begin witli some classical considerations and we focus on 
the gravitational contributions to C, Cj and for C only on the Euclidean piece. For the matter 
contributions and the Lorenzian piece the necessary, completely analogous manipulations can be 
found in [ITj. Consider a partition V of S into cubes □ so that 

^ = J2 f d'a ^f\CP^Q^KJ~C^\{a) (3.7) 

neV 

Let Vo(n) be the coordinate volume of □ in any coordinate system and let cr(n) be some coordinate 
point inside □ with respect to the same coordinate system. Then we can write (13.71) as limit, in 
which the partition becomes the continuum, of the following Riemann sum approximation of the 
above integral 

" = E ^\C'-Q^'C,a\{a{n)) (3.8) 

neP 

Using the classical identities 

where /, J, .. = 1,2,3 label a basis r/ = —iai (where aj denote the Pauli matrices) in su(2) and 
denotes th triad it is not difficult to verify that 

C^ = [Ti{B)]\ Q^^Cfik = [Ti{Br,)f/A=:C] (3.10) 

Here we have introduced the magnetic field B\ = \e^^^Fli and have set = BjTj, ej = CjTj, B = 
B^Cj where ej denotes the cotriad. We may further write 



H = ^im E V\Ciay-S''Cj{n)Cj{n)\ (3.11) 

where 



Des 



C(n) := / rfV C((t), C/(n) := / d'a Ci{a) (3.12) 



The strategy is now to quantise the objects fl3.12p and to define 



p^s 

neS 



H:= lim J2 yl^i^V C{n) - S^JCiiay Cj{a)\ (3.13) 



provided the limit exists. For C(n) this has been done in the literature [TB] and we follow the 
same strategy here. In fact we can treat both C, Cj in a unified way. We have with tq := I2 

f d^aTT{BT^)= [ Tr(FAer^) = i/ Tt{F A {V{n), A} r^) (3.14) 

where 

V{a) = [ £a ^/det{Q) (3.15) 



□ 



is the physical volume of □. Actually there is a sign of det(e) involved in fl3.15p but this is cancelled 
in the squares that appear in (13.111) . 
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The virtue of writing (13.141) in this form is that fl3.15p can be quantised on the LQG Hilbert space, 
hence one replaces the Poisson bracket by the commutator divided by ih. Thus one is left with the 
quantisation of the connection A and its curvature F. This is the source of many ambiguities already 
in unreduced LQG because A, F do not exist as operators, what exists are holonomies along paths 
and loops respectively which can be used in order to approximate A, F respectively. However, while 
classically there are infinitely many ways to do this with the same continuum limit, in the quantum 
theory each choice leads to a different regularised operator in unreduced LQG, see [15]. In unreduced 
LQG one can still argue that most of the uncountably infinite number of choices are gauge related 
under the spatial diffeomorphism group and in fact spatial diffeomorphism invariance is used in order 
to carry out the limit V —>■ X in a specific operator topology [11], |15]. However, in reduced LQG the 
spatial diffeomorphism group is no longer a gauge group, it is a symmetry group of the dynamics. 
Therefore these two arguments are no longer available and therefore the ambiguity issue appears to 
be much worse in reduced LQG. This is the first indication that calls for the AQG generalisation. 

In the next paragraph we will discuss to what extent those ambiguities persist in reduced LQG, 
in the paragraph after that we use the AQG reformulation. 

3.3.2.2 Reduced LQG: Embedded graphs We want to define the Hamiltonian operator H 
on the Gauss invariant Hilbert space of LQG which we will denote by Ti. This Hilbert space has 
an orthonormal basis consisting of spin network functions T^j,/ where 7 is a (semianalytic) graph 
embedded into S, j = {je}eeE{-y) is a collection of non vanishing spin quantum numbers (one for each 
edge) and / = {/j;}j;Gy(7) is a collection of Gauss invariant intertwiners (one for each vertex). There 
is a unitary action of the active diffeomorphisms on this Hilbert space defined by 



In unreduced LQG the diffeomorphisms are considered as gauge transformations and therefore the 
states (13.161) are all gauge related. In the reduced formalism of this paper the states of the form are 
physically distinguishable. Therefore it does not make physical sense to construct diffeomorphism 
invariant distributions which sometimes are used in the construction of Hamiltonian or master con- 
straint operators as already pointed out. 

This last point has crucial bearing on the quantisation strategy: If we want to preserve the 
classical symmetry of the Hamiltonian operator under diffeomorphisms, then this operator must be 
quantised in a graph non changing way [40j on Ti. By this is meant the following: Let T-C^ be the 
closed linear span of spin network states over 7. Then Ti is the direct sum of the T-C^ , that is 



which shows that the physical Hilbert space Ti. is non separable. This is an important difference with 
non reduced LQG where the physical Hilbert space can be made separable if one extends the spatial 
passive diffeomorphism group beyond the differentable category [H] . This is a second indication that 
one should possibly leave the strict realm of (reduced) LQG and pass to another framework where 
non separable Hilbert spaces can be avoided. This calls for the AQG extension [IT] which we discuss 
in the subsequent paragraph. 

In any case, graph non changing in the sense of [IQ] now means that the operator H should 
preserve each TV separately! This appears as if we had to assume an infinite number of conservation 
laws that the classical theory did not have which is a second point to worry about and presents a 
third motivation to switch to the AQG extension of LQG. However, let us see how far we can get 



(3.16) 



(3.17) 
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within the usual formalism. To that end, we use the notion of a minimal loop originally introduced 
in [28| and also used to some extent in 1161. 



Definition 3.1. Given a graph 7, consider a vertex v G V{'~f) and a pair e, e' G -£^(7) of distinct edges 
incident at v and with outgoing orientation. A loop a-y,t,,e,e' in 7 starting at v along e and ending 
at V along (e')""'^ is said to be minimal provided that there exists no other loop in 7 with the same 
properties and fewer edges traversed. The set of minimal loops in 7 with data v, e, e' will be denoted 

by -£7,j;,e,e' • 

Notice that the definition is background independent and diffeomorphism covariant. 
Given a graph 7 and a vertex v G V{j) we define for = 0, 1, 2, 3 

xTi{T^A{a)A{eK)[A{eKy\ (3.18) 

where T„(7) is the set of ordered triples (i.e. order matters) of distinct edges of 7 incident at v taken 
with outgoing orientation, A{p) denotes the holonomy of the connection A along a path p and 



K = 4. / 1 4^ E ^(^1' es) e^MN XM Xg\ (3.19) 

V ei,e2,e3GT„{7) 

is the projection of the volume operator [l2j tc0 TP for an infinitesimal neighbourhood of v. Here 
o'viei, 62, 63) is the sign of the determinant of the matrix formed by the tangents of those three edges 
at V and X^ denotes the right invariant vector field on SU(2) associated with the copy of SU(2) 
coordinatised by A{e). 
Finally we set 



:— E y \P'y[C\,vCy,v — -Cy^^vCi,'y,v]P'y\ (3.20) 

where Pj : Ti 1-C denotes the orthogonal projection and makes sure that H is not graph changing, 
i.e. preserves H."' . The dagger operation is that on Ti for the operator defined in (13.181) using that 
entries of holonomies matrices are just multiplication operators and that V^^^ is self adjoint. 
The operator H is now simply 

H = (3.21) 

It is easy to check that it is diffeomorphism invariant 

f/((^)Hf/((^)-i = H (3.22) 

for all Moreover, it is manifestly Gauss invariant. One may ask what happened to the limit 
V ^ S. The answer is that we define the operator H as in fl3.22p and just check that its expectation 
values with respect to suitable semiclassical states reproduces the classical function H. Such states 
in particular must use sufficiently large and fine graphs in order to fill out S. What the operator 
does on small graphs is irrelevant from the point of view of the classical limit. 

With the methods of ^\ one should be able to verify that on such graphs the semiclassical limit 
of the operator is correct. However, that calculation is of course graph dependent. 



^''The alternative volume operator [13] was ruled out in [33] as inconsistent with the classical Poisson bracket identity 
(|3.14[) . In unreduced LQG one could still say that the volume operator and the Poisson bracket identity are relations 
among non observable objects but this is no longer true in reduced LQG as considerd here and hence the objection 
[33j must be taken seriously. 
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3.3.2.3 Reduced AQG: Abstract (algebraic) graphs One of the motivations for the AQG 
extension of LQG is the graph dependence of the semiclassical calculations. The other is the necessar- 
ily graph perserving feature of diffeomorphism invariant operators which appears to say that there is 
an uncountably infinite number of conservation laws that the classical theory does not have. Finally, 
the non separability of the Hilbert space Ti even if S is compact without boundary is disturbing. 
In a sense, to use all graphs is a vast overcounting of degrees of freedom, at least from the classical 
perspective. To see this, suppose for simplicity that S is topologically (or an open neighbourhood 
thereof) and thus can be covered by a single coordinate system. Consider piecewise analytic paths 
which consist of segments along the coordinate axes. Likewise, consider piecewise analytic surfaces 
which are composed out of segments of coordinate planes. It is clear that the holonomies along those 
kind of paths and fluxes through that kind of surfaces separates the points of the reduced phase 
space. 

It is true that also in canonical QFT the qantum configuration space is always a distributional 
enlargement of the classical configuration space. However, there it is never the case that the label 
set of those fields is uncountable. For instance, in free scalar field theory on Minkowski space the 
quantum configuration space consists of Schwarz distributions rather than smooth functions. The 
label set of the fields consists of test functions of rapid decrease which are dense in the Hilbert space 
of square integrable functions on and there exists a countable orthonormal basis of that Hilbert 
space consisting of Schwarz functions (e.g. Hermite functions times a Gaussian). Thus, the quantum 
fields are tested by a countable set of test functions and an orthonormal basis in the QFT Hilbert 
space is labelled by that countable set. In LQG on the other hand the quantum connections are 
tested by all graphs which is an uncountable set and states over different graphs are orthogonal. So 
the situation is completely different which seems to be the price of having a diffeomorphism covariant 
theory [nHH]. 

One could of course restrict the labels to those mentioned above but these would not be preserved 
by diffeomorphisms. It is true that the diffeomorphic image of a coordinate segment can be approxi- 
mated by coordinate segments, however, the length of say a rotated segement when approximated by 
a staircase will differ largely from the original length. The same happens for areas of surfaces. The 
only chance that this does not happen is for observables that are integrals over three dimensional 
regions as pointed out in [45] . 

To make progress on those issues we therefore will restrict attention to operators that come from 
integrals over regions of S such as the volume operator or the Hamiltonian operator. This does not 
mean that one cannot construct length and volume operators, one just has to define them in an 
indirect way, see [16]. In fact, we will only consider quantising functions which are Diff(iS) invariant. 
The motivation for doing this is that in physics we do not specify spatial regions by considering a 3D 
subset R oi S and define, say, a Diff(A') invariant volume functions (i.e. a function invariant under 
passive diffeomorphism starting from the unreduced formalism) by 



where xr denotes the characteristic function of the set R. Rather we use observable matter for doing 
this. To be sure, f l3.23p is Diff(A') invariant, being the integral of a scalar density over all of X. In 
fact, we can pull back this expression to S and obtain 



where q = {S ^) * q which would be a mathematically natural object to consider in the reduced 
theory (after further applying a^o(.)^o=T-_7-). It is, however, not DifT(5) invariant. However, from 




(3.23) 



V{R) = [ (fa Xi?(fT) v/det(g)(a) = / (fa ^dei{q){a) 



(3.24) 
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the point of view of observation one would rather hke to consider an object of the form 



Vil) := / d^x xMx)) ^/detiq){x) (3.25) 



where / is a subset of the real axis and is a scalar field. Notice that f l3.25p is Diff(A') invariant but 
it is not a Dirac onbservable yet. It measures the volume of the subset of X in which has range in 
I. Now we apply the map O'^ and obtain immediately 

Ov(/)(r,a) := ^ rfVx/($(r,a)) ^det (Q)(r, a) (3.26) 

where $(r, a) = 0^x){'t,o') (for any x) is the Dirac observable associated to 0. Curiously, f l3.26p is 
a Dirac observable and it is Diff(i5) invariant. It measures the physical volume of the region in S 
where the physical scalar field $ ranges in /. The argument shows that Diff(iS) invariant observables 
naturally arise from the point of the unreduced theory and from operational considerations. 

Having motivated to consider only Diff(iS) invariant observables we are now ready to consider the 
AQG framework. Since for such observables the coordinate system plays no role we generalise from 
embedded to non embdeed graphs and the above argument shows that infinite cubic algebraic graphs 
should be sufficient although a generalisation to arbitary countable algebraic graphs as sketched in 
[17] would be desirablj^. In this paper we will just consider the cubic graph for simplicity. 

At the algebraic level the notion of Diff(iS) and even of S itself is meaningless. Notice that in AQG 
the infinite algebraic graph is a fundamental object. This fundamental graph does noch change. What 
does change under the dynamics are subgraphs of the algebraic graph. In other words, subgraphs of 
the fundamental algebraic graph are not preserved under the quantum dynamic^. The definition of 
H in AQG is much simpler and no longer involves the projection operators P^, so we do not have the 
awkward conservation laws any longer. In fact, there is no dependence on any algebraic subgraph 
whatsoever. In complete analog}0 to p/T] it is given by the following list of formulae 



^ Si,S2,S3=±l 

xTT{T^A{a^.I^suhs2)M(^v;l3S3)[M'^v;l3S3)~^,Vv]) (3.27) 



where e^-js is the edge beginning at v in positive (s = 1) or negative (s = —1) / direction and ay.is,js' 
is the unique minimal loop in the cubic algebraic graph with data v, ev-is,^v;js'- Formula fl3.27p is 
actually the specialisation of f l3.18p . The operator is the algebraic volume operator 




^ V siS2S:,eiJ^ eLMN XL ^ Xf^ | (3.28) 

AO / J ^v\1sy ^';js2 *^v\Ks3 i ^ ^ 

S\,S2,S3=±\ 



^'"^The idea would be to consider the most general such graph which is the maximal algebraic graph. This is an 
algebraic graph with a countably infinite number of vertices and with a countably infinite number of edges between 
each pair of vertices including loops. This generalises the notion of a complete graph which is a graph in which a 
single edge connects each pair of vertices. 

^^This bears some resemblance with the models for emergent gravity considered in |47j although the dynamics of 
those models not obviously models the dynamics of H. 

^^In we considered the extended Master constraint which involves, in the language of this paper, [C^ + 
Q^^GjCk\|^JAQi[ff) rather than ^\C'^ - QJ'=CjCfe|. Hence apart fr om the sign in front of Q^^CjCk we only need 
to change the power of the volume operator from Vv in |17j to Vv here. 
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Finally 




(3.29) 



V 



where the sum is over all of the infinite number of vertices of the algebraic graph. The operator 
f l3.29p is manifestly Gauss invariant. 

The Hilbert space of AQG is the infinite tensor product (ITP) of Hilbert spaces L2{SU{2), dfin), 
one for each edge of the graph (this can be generalised to defining different ITP's that come into 
play when constructing Gauss invariant states). This Hilbert space is not separable but it is a direct 
sum of separable Hilbert spaces which assume a Fock like structure and which are preservecQ by 
H. As far as the symmetry group is concerned, at the algebraic level for instance we no longer 
have spatial diffeomorphisms. However, we have its algebraic version which consists in the following: 
Consider the master constraint like functional 



where 6 > a > are any real numbers. Then a classical function F is invariant with respect to 
transformations generated by H,Cj respectively if and only if {F, {F, M}}m=o = 0. The functional 
(I3.30p can be quantised on the AQG Hilbert space by literally the same techniques as in [TT]. Thus, 
we have the possibihty to analyse the anomaly issue with respect to (S at the algebraic level as well. 

Finally, we can compute the expectation value of H with respect to semiclassical states as in [17] 
and to zeroth order in h we should find that the classical value is reproduced with small fiuctua- 
tions. As for the master constraint, the really astonishing fact is that H is a finite operator without 
renormalisation thanks to our manifestly background independent formulation. Namely, at the fun- 
damental quantum level the operator algebra is labelled by a single, countably infinite abstract, that 
is non embedded, graph F. There is no such thing as a lattice distance which would need a back- 
ground metric. However, the semiclassical states depend on a differential manifold X, an embedding 
Y of the algebraic graph F into X, a cell complex y (F)* dual to Y{T) as well as a point {Aq, Eq) in 
the classical reduced phase space. Thus, the semiclassical states make contact to the usual (reduced) 
LQG formulation which in particular uses an at least topological manifold X. Hence AQG describes 
all topologies simultaneously. The point is now that, since F is an infinite graph, the embedding of 
F can be as fine as we wish, with respect to the spatial geometry described by Eo even if X is not 
compact^. The expectation values of our operators such as H will now give, to zeroth order in h, a 
Riemann sum approximation of the desired continuum integral H as in fl3.14p in terms of holonomies 
along edges of the embedded graph and the volume of the cubes in the dual cell complex. That 
Riemann sum will approximate the integral the better, the finer the embedding. It is in this sense 
that in (13.291) no continuum limit has to be performed. 



^^We do not know whether H is densely defined on all of the ITP. However, if it is defined on a single vector in a 
given separable sector then it is densely defined on the entire sector. Now each separable sector of the ITP is labelled 
by a cyclic vector which is explicitly known. Now H is defined on a given H. if and only if it is densely defined on 
the corresponding sector. Hence, for each we just have to perform this test and we simply remove the sectors from 
the ITP on which H it is not densely defined, if any, since they are unphysical. H is certainly densely defined on the 
sectors built from semiclassical ft, hence the surviving part of the ITP certainly includes all the semiclassical states. 

^^In the compact case the embedding necessarily has accumulation points but we can choose our states not to be 
excited on edges that are mapped under the embedding into a suitably small neighbourhood of every accumulation 
point. 




(3.30) 
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4 Summary and Outlook 



As compared to the Master Constraint Programme [16] the present framework has the advantage 
that the Master constraint and its solutions are not needed. We directly consider a representation 
of the gauge invariant phase space and its Hamiltonian. Celebrated results of unreduced LQG such 
as the discreteness of the spectrum of kinematical geometric operators [l2l HSl HH |19] which is not 
granted to survive when passing to the physical Hilbert space [501 [5T] in the usual Dirac constraint 
quantisation now becomes a physcial prediction if the curves, surfaces and regions that one measures 
length, area and volume of are labelled by dust space. The Gauss invariant [52] kinematical coherent 
states [18] of unreduced LQG now become physical coherent states. 

However, the physical Hilbert space of reduced LQG is non separable which appears to be a vast 
overcounting of quantum degrees of freedom. Passing to AQG means to switch from embedded to non 
embedded graphs and thus removes the overcounting. Since for spatially diffeomorphism invariant 
operators (on dust space) such as the Hamiltonian H or any other operationally interesting observable 
(which does not refer directly to the dust label space) the embedding of a graph is immaterial, we 
can consider the AQG reformulation as an economic description of reduced LQG in the sense that 
diffeomorphism related embeddings would lead to isomorphic sectors superselected by this kind of 
observables. The additional advantage of AQG is that it does not require a topological manifold and 
that it is free from complications that have to do with graph preservation. 

The challenge of the present framework is to implement the (algebraic version of) the symmetry 
group & in the definition of H which will require tools from lattice gauge theory. The final AQG 
version of the reduced phase space is in any case very similar to Hamiltonian lattice gauge theory 
with the important difference that no continuum limit has to be taken which is why the theory is 
UV finite. Another important question is how one can understand from the complicated, non per- 
turbative Hamiltonian H the significance of the standard model Hamiltonian on Minkowski space. 
The answer to this question must lie in the construction of a minimum energy eigenstate of H which 
is simultaneously a minimal uncertainy state Q for all the observables and which is peaked around 
flat vacuum (no excitations of observable matter) spacetime. Presumably, if one studies matter ex- 
citations of Q and considers matrix elements of H in such states then the resulting matrix elements 
can be considered as the matrix elements of an effective matter Hamiltonian on Minkowski space 
which should be close the Hamiltonian of the standard model on Minkowski space. This expectation 
is supported by the analysis of [6l [7] which shows that the equations of motion of the gauge invariant 
geometry and matter degrees of freedom perturbed around a homogeneous and isotropic (FRW) 
solution is described effectively by the usual Hamiltonian on a FRW background expanded to second 
order in the perturbations. Of course, this is only a classical argument. See [2B] for more details 
about the quantum aspects of this idea. We leave this and the research projects mentioned in the 
introduction for future analysis. 
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